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Abstract 

We investigate the distribution of the time spent by a random walker to the 
right of a boundary moving with constant velocity v. For the continuous-time 
problem (Brownian motion), we provide a simple alternative proof of Newman's 
recent result [J. Phys. A 34, L89 (2001)] using a method due to Kac. We then 
discuss the same problem for the case of a random walk in discrete time with an 
arbitrary distribution of steps, taking advantage of the general set of results due 
to Sparre Andersen. For the binomial random walk we analyse the corrections to 
the continuum limit on the example of the mean occupation time. The case of 
Cauchy-distributed steps is also studied. 



To be submitted for publication to Journal of Physics A 
P.A.C.S.: 02.50.Ey, 02.50. Ga, 05.40.+j 

1 godreche@spec. saclay.cea.fr 
2 luck@spht. saclay.cea.fr 
3 URA 2306 of CNRS 



1 



1 Introduction 



Consider a Brownian particle, starting from the origin, whose position x t satisfies the 
Langevin equation 

dx t 

where rj t is Gaussian white noise, such that (rj t ) = and (f]tVt') — 2D 5(t — £'). 

In a recent paper Newman considered the following question: What is the distri- 
bution of the length of time spent by the particle to the right of a boundary moving with 
constant velocity vl This quantity, denoted by T 4 + , is known as the occupation time of 
the half-line located to the right of the boundary. It reads 

ft 



t+= f dt'i t ,, i t , = e(x t ,-vt'), 

Jo 



where Q(x) is the Heaviside function. The indicator random variable It is therefore 
equal to 1 if x t > vt, and else. Similarly, the occupation time to the left of the moving 
boundary is denoted by Tf, and such that T t + + Tf = t. 

A number of past studies as well as more recent ones have been devoted to the 
statistics of the occupation time of stochastic processes, either in probability theory [Q, 
§, § §, or in statistical physics §, 0, §, |, 0, [TT], |TJ, |TJ, |T5| in the context of 
persistence. 

A derivation of the probability density of T t + , f T +(t,r) = (YP{T^ < r)/dr, is given 
in Ref. U, with the result 

f Tt+ (t,r)=F + (r,v)F-(t-T,v), (1) 

where 



and erfc is the complementary error function. 

In particular, in the case of a static boundary (v — 0), we have 

/ T+ (t,r) = - / (3) 

hence the fraction of time T t + /t spent by the Brownian particle to the right of the origin 
admits a limiting distribution as t — > oo, which reads 

lim f T +Jx) = — 1 (0 < x < 1). (4) 



t— >oo 



7T\ X[l — X) 



The arcsine law |l^| is thus recovered. 

The aim of this note is to complement Newman's work in two directions. 
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Firstly, we give an alternative, simpler derivation of equations ([!]), @, using a method 
originally due to Kac ||. 

Secondly, we discuss the corresponding problem for a sum of random variables, i.e., for 
a random walk in discrete time, with an arbitrary distribution of steps, either narrow or 
broad. The general set of results due to Sparre Andersen |Tj], [RJ is the starting point of 
this analysis. A factorisation property of the distribution of the occupation time similar 
to (|I|) holds, the role of F + (t,v) being played by the quantity F^(v), which is simply 
the survival probability of the walk in the presence of the boundary. We investigate two 
examples in more detail. For the binomial random walk, for which a detailed study of 
F^(v) can be found in we analyse the corrections to the continuum limit on the 
example of the mean occupation time; for a Cauchy distribution of steps, we determine 
the probability distribution of the occupation time. 



2 Brownian motion 



The problem of a particle executing symmetric Brownian motion in the presence of a 
boundary moving with constant velocity v is equivalent to that of biased Brownian motion 
with velocity — v in the presence of a fixed boundary, located at the origin. Consider the 
probability for this biased Brownian walk to be at position x at time t, and to have spent 
a length of time equal to r to the right of the origin. The joint probability density of the 
event (x t = x, T t + = r) is denoted by p(t, r, x). We then have 



/oo 
dxp(t,r, 
-oo 



X 



(5) 



The method of KacQ consists in writing a master equation for p(t, r, x). It is an easy 
matter to realize that the latter reads, in the present case, 



dt +U[X) dr 



^d 2 p dp 

dx 2 dx ' 



(6) 



with initial condition p(0, r, x ) = 6(x)5(t). In Laplace space, setting 



p(s,u,x) =£ p(t,r,x), 

t,T 



equation @ yields 



(s + uQ(x))p- D 



d 2 p dp , . 



(7) 



dx 2 dx 

This inhomogeneous differential equation is easily solved. By requiring p(s, u, x) to vanish 
at infinity (x — ► ±oo), we have 

x 

2D. 



p(s, u, x) = A(s, u) x 




Wis 



(-V + y/v 2 + 4Ds 



U , 



(x > 0), 



X 



2D 



(x < 0). 



lr This method has recently been applied to the investigation of the distribution of the occupation 
time of subordinated Brownian motion [Q3| (see also fl4|). 
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The amplitude A(s,u) is determined by the right side of (0), yielding the condition 

dp(t,r,0+) _ dp(t,T,0~) 1 
dx dx D ' 

hence 



v 2 + 4D(s + it) - V^ 2 + 4£>s 
v ' ; 2Du 
Finally equation fl5|) yields 

f T +(s,u) =C (e~ uTt ) = / dxp(s,u,x) = + u, w), (8) 

with 

F (u, v) = - . 9 

These functions are precisely the Laplace transforms with respect to r of -F ± (r, v) given 
in equation (|2[). An inverse Laplace transformation finally leads to the result ([!]). 

Let us now comment on the form of the solution (0). A first comment concerns 
symmetry. It is clear that, for a given velocity v, T^{v) has the same distribution as 
Tf(—v). Therefore 

f T + (_„) (*, T ) = /t- ( B ) (* > r ) = U+ (v) (*> * - r ) • 

This requirement is satisfied by the solution ([!]), since -F ± (t, — 1>) = F^(r,v), according 
to (g), or equivalently ^(u, —v) = F T (u,v), according to A second comment is on 
normalisation. Integrating equation ([!]) upon r G [0, t) yields 

1 = V(T+ <t)= /"dryWfcr) = * F" (*,«), 

the star in the right side denoting a convolution product. Hence in Laplace space the 
equality F + (s, v )F~(s, v ) = 1/s should hold. This is indeed the case, as can be seen 
from (|9|). 

For v = 0, the distribution of the occupation time is given by (Q). The first correction 
to this behaviour for small v is given by 

, , , 1 1 v ( 1 1 \ 

f T +(t, T = + — _ — - - + • • • 

As a consequence, 



showing that the presence of any bias v ^ is relevant in the long-time regime. There is 
actually a non-trivial limiting distribution for T 4 + as t — > oo if v > (and, by symmetry, 
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for T t if v < 0). We have indeed lim^oo F (t,v) = lim s _^ S F ( s , v ) — v/\fD, so that, 
by ®, 

f T+ (r) = ^F + (r,v), (10) 



with the notation T + = lim^oo T 4 + . Using the asymptotic expansion 

r , x e_a;2 / 1 
erfc(x) = — ^ 1 - — + • • 

we see that the distribution (|To"|) falls off exponentially for large r, as 



M(r) ^V^ exp H^)' 

so that all the moments ((T + ) n ) are finite. The latter can be computed by noting that 

fr+(u) = -J=F + (u,v), 

which, expanded around u — 0, leads to 

3 Random walk in discrete time 

Consider a discrete random walk defined by a sum of independent, identically distributed 
random variables: 

n 
i=l 

with an arbitrary distribution of the steps rji (either discrete or continuous, narrow or 
broad). For this random walk, the occupation time to the right of the boundary moving 
with velocity v is defined as 

n 

T+ = ^ Imi Ira = ©(^m ~ V17l) , 

m=l 

hence the indicator random variable I m = 1 if x m > vm, or else. The occupation time 
T~ to the left of the boundary is defined likewise, and such that T+ + T~ = n. As above, 
we note the equivalence of the problem thus stated with that of the occupation time of 
a random walk with biased steps rji — v in the presence of a fixed boundary, located at 
the origin. 



For sums of independent random variables, a result due to Sparre Andersen |17|, |16 
expresses the probability distribution of T+ as the product 

V(T+ = k)= V(T+ = k)V(T~_ k = n-k). (12) 
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In this equation, V(T£ = k), hereafter denoted by F£(v), is the probability that the 
walk remained to the right of the boundary up to time k, or 



F+(v)=V(T+ = k) = (hI 2 ...I k ). 

Similarly, 

F^v) = V(T k — k) — ((1 - A) (1 - h) ... (1 - h)) 

is the probability that the random walker remained to the left of the boundary up to 
time k. In other words, the quantities F^{v) are survival probabilities of the walk in the 
presence of the boundary, up to time k || [7], ^, [18], |19| . For instance, 

Fk ( v ) = 'P ( x m > vm for 1 < m < k) . 

The generating function of the F£(v) is related to the generating function of the 
one-time probabilities (J n ) = V(x n > vn) by |17|, [16 



oo n 



£F+(^ = exp £-</„> • (13) 



n=0 \n=l U 



Equation (0) is the discrete counterpart of equation (|l|). Together with flHf), it provides 
the answer to the question posed (in terms of the one-time quantities (J ra )). 

We now illustrate the above formalism by two examples. First, for a narrow distribu- 
tion of steps 7]i, we determine the continuum limit of equations (|T2|), ([131), thus recovering 
the results (0), (0) obtained for Brownian motion. For the binomial random walk, we 
investigate the corrections to the continuum limit on the example of the mean occupa- 
tion time. Then, for a Cauchy distribution of steps, we determine the distribution of the 
occupation time V{T+ = k) explicitly. 

Consider a narrow distribution of steps, with (rj) = and (rj 2 ) = a 2 . The continuum 
limit is defined as n — > oo, v — > 0, with v-Jnjo = £ fixed. The central limit theorem 
yields 



(4) = V(x n >vn)=V f^= due- 2 / 2 = Ul 



erf^y. (14) 



Let us analyse equation (|T^) in the same limit. Setting z = e~ s , to leading order 
as s — > 0, identifying generating series with Laplace transforms, we obtain the following 
estimates: 

oo 
n=0 



n=1 n ay/2 Jo t ay/2 
Hence finally 




F + {s,v) = -j=====——. (15) 
\/v z + 2a z u + v 
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With the identification a 2 = 2D, the result (||) is recovered. 

For the binomial random walk, the survival probability F^{v) is a highly non-trivial 
function of the velocity v , depending on whether v is rational or irrational, because of the 
underlying lattice structure [|T8|| . In particular, the limit survival probability F + (v) = 
lim^oo F+(t>), which is non-zero for v < 0, is discontinuous at any rational value of 
v. The corresponding discontinuities are algebraic numbers which can be determined 
explicitly fllS] . In the continuum limit (v — > _ ), F + (v) is simply given by \v\y/2/a. 

In order to better understand the nature of the corrections to the continuum limit, 
we consider the simple case of the asymptotic mean occupation time (T + ). For the 
symmetric binomial random walk, and for fixed v > 0, we have 

oo oo n /n\ 

(T + ) = Y,(in) = Y,z- n E (?), as) 

n=l n=l k=k (v) W 

with ko(v) = Int(n(l + v)/2) + 1, where Int(x) is the integer part of x, i.e., the largest 
integer less than or equal to x. As shown in the appendix, the behaviour of this expression 
as v — > is given by 

with 




A = 4 - Cl — = -0.165869209. 



The first term in (|T7|) corresponds to the continuum-limit result (see equation (|TT|) ) 

1 



I Brown 2,1) 

because D = 1/2 for the binomial random walk. The second term in (\Tf\j is surprising 
in several respects: it is of relative order t> 3//2 , instead of the naturally expected v 2 , and 
the corresponding amplitude A is transcendental. 

Figure 1 shows a plot of the difference between the exact value of (T + ), obtained by 
evaluating numerically ([16]), and its continuum-limit expression (fL8|). The sum of the 



last two terms in (17); shown as a thick line, correctly describes the mean asymptotic 
behaviour of the plotted quantity. Superimposed periodic oscillations, with period two, 
are also clearly visible. Similar oscillations, due to the lattice underlying the discrete walk, 
are also encountered when considering other quantities [0 (see especially Figures 14 
and 15 therein). 

Consider finally the case where the steps have a Cauchy distribution 

1 

p(v) 



7r(l + r/ 2 )' 

Because of the stability of the Cauchy law, the probability (J n ) = V(x n > vn) is inde- 
pendent of n: 

r°° 11 

(/„) = V{x n > vn) = a(v) = / dup(u) = arctant>. (19) 

Jv 2 ix 
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C\2 

> 



V 




Figure 1: Plot of the difference between the mean occupation time (T + ) for the symmetric 
binomial random walk, and its continuum- limit expression (|lg|), against 1/v. Thick line: sum 
of last two terms of equation ([l7]) . 



Hereafter we denote this expression by a, for short. Using (|13D, we get 



n=0 



hence 



with 



F+(v) 



T(n + a) _ n~ e ^ 
n\T(a) ~ T(a) 



(n > r 



6{v) = 1 — a = — I — arctan-y. 

2 7T 



(20) 



The survival probability F+(v) falls off as a power law, with a continuous family of 
persistence exponents 6{y ) |T|, |5|, |K) . 



From equation (0), we obtain 

T(a + /c) r(n + 1 — A; — a) (n\B(a + k,l — a + n — k) 



v(T+ = k 



k\T(a) (n-k)\T(l-a) 



B fa, 1 — a) 



(21) 
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where the beta function is defined as 



flM ) = /W>(i-«)-' = ^J^ 

Jo T (a + b) 

Defining the (3 density as 



Mx) = 1 ^- 1 TX a - 1 (l - x) b - 1 (0 < X < 1), 
B{a, b) 

one can rewrite (2~I) as 

7>(T+ = fc) = J\ul^u k {l-u) n ~ k (3 a ^ a {u). (22) 
In the continuum limit where n and are simultaneously large, with a fixed ratio 



x = k/n, the binomial distribution inside (22) converges to 5(u — x), so that the limiting 
probability density function of T+ fn reads 



lim f T+/n (x) = (3 a ^ a {x) = _ J)-. ( 23 ) 



In particular, if t> = 0, then a = 1/2, and one recovers the arcsine law (§). 

To conclude, let us briefly consider the case of a moving boundary whose position 



obeys an arbitrary power law: X(t) =wt v |2(| |2l], [ll 

If the steps have a narrow distribution, the continuum description can again be used 
in the regime of long times it ^> 1) and weak bias (|w| «C 1). A generalisation of flU] ) 
shows that the mean occupation time reads 

The case v = 1/2 therefore demarcates between two regimes. For v < 1/2, the bias w is 
irrelevant, so that (T t + ^ ~ t/2, and the arcsine law Q) still holds asymptotically. To the 
contrary, for v > 1/2, any weak bias is relevant, so that a non-trivial limiting law for the 
occupation time T + is expected (for w > 0), generalising (|10|), with (T + ) ~ (a /w) 2 ^ 2u ^ 1 \ 
In the marginal situation of a parabolic boundary {y = 1/2), there is a continuously 
varying persistence exponent 9(w) |20|, ^TJ, and the fraction T + /t admits a non-trivial 
limit distribution, which also continuously depends on w. 

If the steps have a symmetric broad (Levy) distribution, with tails falling off as 
p{rf) ~ > with < fi < 2, the above discussion on the relevance of the bias w 

still applies, with the marginal situation being v — \j \i. The case of Cauchy-distributed 
steps in the presence of a ballistic boundary (// = v = 1) is an interesting example 
of this marginal situation, where the dependence on the bias w = v of the persistence 
exponent (|20|) and of the limit distribution (B3|) are known explicitly. 
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Appendix: Expansion for v — > of expression 



In this appendix we investigate the behaviour as v — > of expression (|T^) of the mean 
occupation time (T + ), i.e., 

co co n / \ 

n=l n=l k=k (v) V / 

We set m = 1/v, and introduce the Laplace transform L(s) =£ (T + ). The above expres- 

u 

sion yields 

with fc = fc o (0) = Int(ra/2) + 1. 

Introducing the contour-integral representation 

(n\ r dz 1 



\k) J 2m z k+1 {l - z ) n - k+v 

where the contour encircles the origin, and summing over k at fixed £ = 2k — n > 1, we 
obtain 

L(s ) = ^ll^y i 

s J 2vri ^ (2zy(Az(l - z) - e" 2 ^) ' 

hence, after some algebra, 



s J^i w i s fr{ Wz 



with W e = VI - e~ 2s/e - 

In the regime of interest (s — > 0), we have 

L(s) = L (s) + L^s) + ■ ■ ■ , 

with 

1 °° CO / 1 1 ' 

L (s) = V Vie-^, L x {s) = V - - 

and so on. 

The leading series L Q (s) has to be investigated in some detail, by means of its Mellin 
transform M (x), for which we obtain a closed- form expression: 

M (x) = dss x - 1 L (s) = 2 2 " x r(2x-3)C(a;-2) (Rex > 3), 
J o 

where ( is Riemann's zeta function. The behaviour of the subleading series Li(s) can be 
estimated to leading order, replacing the sum over £ by an integral over V2s£. We thus 
obtain Li(s) « l/(3s). 
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Inverting successively the Mellin and Laplace transforms, we obtain expression 

i.e., 

T + ) = 1 1 1 

with 

A = J- ((--) = -0.165869209, 



71 V 2, 

and where the dots stand for a contribution going to zero as v — > 0. 
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